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LONG  TERM  GOALS 

Time  and  space  dependent  radiance  distributions  at  the  sea  surface  are  a  function  of  the  shape  of  the 
incident  distribution  on  the  surface,  modifieation  by  the  sea  surface  itself  from  topography  and 
transmission  charaeteristies,  and  alteration  by  the  Inherent  Optical  Properties  (lOPs)  of  the  surfaee 
ocean.  The  long  term  goal  of  the  proposed  work  is  understanding  this  last  controlling  factor.  With  a 
knowledge  of  the  lOPs,  radiance  fields  ean  be  directly  computed  from  the  ineident  field  using  the 
equation  of  radiative  transfer,  now  embedded  in  commereially  available  code  (e.g.,  Hydrolight). 

With  the  state  of  eurrent  teehnology  and  methodologies,  the  primary  obstacles  in  understanding 
subsurface  lOPs  and  their  high-frequeney  dynamies  are  a  laek  of  1)  volume  seattering  instrumentation, 
2)  comprehensive  inversion  models  linking  the  lOPs  with  the  causative  bubble,  partieulate,  and 
dissolved  matter  in  the  water  (whieh  in  many  eases  will  require  input  dependent  on  1),  and  3)  suitably 
stable,  non-intrusive  platforms  to  sample  the  subsurface  oeean.  The  first  two  challenges  are  addressed 
in  this  projeet. 

OBJECTIVE 

There  are  two  overarching  objectives  for  this  projeet: 

1)  To  develop  an  in-situ  volume  seattering  function  device  measuring  volume  seattering  from  10°  to 
170°  at  10°  intervals  and  sampling  rates  of  1  s'^  or  better  to  sample  the  VSF  in  near-surface  waters;  and 

2)  To  develop  and  refine  lOP  inversion  models  to  resolve  partiele  field  charaeteristies  on  small  spatial 
(cm’s)  and  temporal  (<1  s)  seales  in  near-surfaee  waters. 
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APPROACH 


Our  design  for  the  VSF  deviee  is  illustrated  in  Figure  1.  The  deviee  is  ealled  MASCOT  (Multi- Angle 
scattering  Optieal  Tool).  The  souree  beam  is  a  30  mW  650  nm  laser  diode  expanded  with  a  Gallilean 
2X  beam  expander  to  an  approximately  3mm  X  8  mm  elliptieal  shape.  A  wedge  depolarizer  is  used  to 
provide  the  unpolarized  light  needed  for  VSF  determinations.  Seventeen  independent  silieon  diode 
deteetors  spaeed  in  a  semieirele  10  em  around  the  sample  volume  measure  volume  seattering  from  10° 
to  170°  at  10°  intervals.  The  total  pathlength  for  all  seattering  measurements  is  20  em.  Independent 
deteetors  allow  resolution  of  the  VSF  without  any  moving  parts  and  time-eonsuming  seanning. 
Additionally,  eaeh  deteetor  ean  be  optimized  for  its  speeifie  measurement.  Deteetor  field-of-views 
(FOVs)  range  from  0.8°  to  5°  for  the  different  deteetors,  with  the  narrowest  FOVs  assoeiated  with  the 
deteetors  measuring  seattering  in  the  forward  direetion.  Using  proprietary  eleetronies,  a  20  Hz 
sampling  rate  for  all  ehannels  has  been  aehieved  while  maintaining  a  worst  ease  signahnoise  of  300:1. 
Relatively  fast  sampling  rates  are  important  in  resolving  VSFs  in  the  highly  dynamie  oeean  subsurfaee. 

MASCOT  ealibration  parameters  have  been  derived  in  solutions  of  mierospherieal  beads  and  the 
partiele  standard  Arizona  Road  Dust  (AZRD;  Powder  Teeh.  Ine.).  A  referenee  deteetor  may  be  used  to 
aeeount  for  any  fluetuations  in  souree  intensity.  Coineident  measurements  of  beam  attenuation  are 
required  to  eorreet  for  light  losses  along  the  optieal  path  of  the  seattering  measurements  (all 
pathlengths  nominally  -'20  em).  Algorithms  have  also  been  eontemplated  that  would  aeeount  for 
attenuation  without  a  separate  direet  measurement.  The  design  in  Fig.  1  is  made  as  small  as  possible 
while  aeeommodating  all  neeessary  hardware,  a  eriterion  established  in  part  to  minimize  this  effeet. 
Two  MASCOT  prototypes  are  being  developed  for  redundaney  in  ease  a  problem  with  a  eomponent 
arises.  Two  systems  allows  for  eoneurrent  testing  of  different  sensor  eomponents. 

For  the  inversion  modeling,  we  will  extend  the  eapabilities  of  existing  models  (Twardowski  et  al. 

2001;  Twardowski  and  Zaneveld,  2004)  by  ineorporating  input  from  new  VSF  measurements  and  by 
adding  bubble  partiele  populations  (elean  and  eoated)  in  the  models.  At  a  minimum,  we  expeet  the 
inversion  model  will  give  us  eomponent  eoneentrations  and  size  distributions,  with  the  eomponents 
being  organie  partieles  (living+detrital),  minerals,  and  bubbles.  We  expeet  to  eo-deploy  the  MASCOT 
and  the  eommereially  available  near-forward  VSF  deviee  LISST  (Sequoia  Ine.)  in  order  to  eapture  the 
VSF  with  good  resolution  from  -0.1  degrees  to  170  degrees. 

RaDyO  field  deployments  are  eurrently  planned  in  eollaboration  with  other  RaDyO  investigators  off 
Seripps  Pier  in  January  2008,  in  Santa  Barbara  Channel  in  September  2008,  and  off  the  Hawaii  islands 
in  May  2009.  For  these  efforts,  we  plan  to  eoneurrently  deploy  CTD  +  AC9  (+  ACS)  +  MASCOT  + 
LISST  +  ECO  sensors  in  a  real-time  vertieal  profiling  system  (Fig,  2).  We  also  expeet  to  integrate  two 
additional  sensors  for  these  deployments:  1)  a  high  sampling  rate  fish-eye  lens  radiometer  eurrently 
being  developed  by  Marlon  Lewis  and  Seott  MeLean  of  Satlantie,  and  2)  a  bubble  aeoustie  resonator 
provided  by  David  Farmer  and  Svein  Vagle.  For  the  Seripps  Pier  experiment,  a  mobile  deployment 
platform  will  be  eonfigured  so  that  measurements  may  be  made  along  the  length  of  the  pier,  with  the 
expeetation  that  higher  eoneentrations  of  bubbles  will  be  observed  nearing  the  surf  zone.  For 
subsequent  field  work  in  SBC  and  Hawaii,  we  will  deploy  our  sensor  paekage  off  the  RW  Kilo  Moana. 
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Figure  1.  Oblique  view  illustration  of  the  MASCOT.  The  VSF  is  resolved  from  10  to 
1 70  degrees  in  10  degree  intervals.  Detectors  are  wedge  shaped  and  arranged  in  a  semi¬ 
circle  on  an  aluminum  frame  to  minimize  reflections  and  perturbation  of  the  water 
sample  in  the  remote  volume  (center  of  semi-circle).  The  source  assembly  includes  a 
30  mW 650  nm  laser  diode,  reference  detector,  beam  expander,  and  wedge  depolarizer. 
Wiring  from  all  the  detector  modules  and  the  source  module  feeds 
to  a  data  handling  unit. 


WORK  COMPLETED 

•  The  first  MASCOT  prototype  was  completed  in  October  2006; 

•  A  summary  of  calibration,  field  deployment,  and  data  processing  activity  is  provided  in  Table  1 
below; 

Table  1.  Timeline  of  MASCOT  activity  in  FY07.  Time  after  May  2007 
has  been  devoted  to  calibration  analyses  and  data  processing. 


OCTOBER 

week  of  the  9th:  delivery  to  WLE 

18-22:  Deployment  in  Pen  River,  ME  (observed  bubbles) 

DECEMBER 

14-15:  Calibration  exercise;  MASCOT  and  an  AC9  in  series  of  suspensions  of 
Arizona  Road  Dust  (AZRD)  and  microspherical  beads  of  1 .992+/-0.025 
um  size. 

JANUARY 

rework  at  WL  HQ:  dark  count  offsets  between  air  -  water  improved 

FEBRUARY 

6-7:  Calibration  exercise;  MASCOT  and  an  AC9  in  series  of  suspensions  of 

Arizona  Road  Dust  (AZRD)  and  microspherical  beads  of  1 .992+/-0.025 
um  size;  this  time  with  all  WLE  ECO  BB  sensors  as  well. 

MARCH 

First  2  weeks:  Deployment  in  Hawaii  (SORTIE) 

21-28:  Deployment  in  northwest  Atlantic  (OCVAL:  Endeavor) 

MAY 

4-10:  Deployment  in  Hudson  plume  region  (OCVAL) 
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RESULTS 


Key  results  in  FY07  were  1)  the  MASCOT  prototype  was  fabrieated,  tuned  and  tested;  2)  a  viable 
ealibration  algorithm  for  the  MASCOT  was  developed;  and  3)  a  substantial  body  of  field  data  was 
eolleeted  from  several  diverse  water  types  around  the  world.  Fig,  2  shows  the  MASCOT  deviee  in  a 
eustom  eage  with  additional  optieal  sensors  ready  for  deployment.  Below  is  a  deseription  of  the 
calibration  method  that  was  developed  and  some  preliminary  field  results. 

To  calibrate  each  of  the  17  raw  signals  to  a  volume  scattering  coefficient,  P(0)  (m'^  sf'),  a  dark  offset 
(DO)  must  be  subtracted  from  the  raw  digital  counts  and  the  result  multiplied  by  a  scaling  factor  (SF). 
Dark  offset  variability  improved  after  the  January  2007  rework.  The  change  in  dark  offset  counts 
collected  in  air  versus  water  ranged  between  10%  and  310%  of  the  approximate  count  levels  measured 
in  the  lab  in  a  c=0.1  m'^  solution  of  the  particle  standard  Arizona  Road  Dust  (AZRD).  A  good  target  is 
a  few  percent  for  all  channels.  After  the  rework,  the  change  in  dark  offset  counts  between  air  and 
freshwater  ranged  between  6%  and  160%  of  these  count  levels.  For  all  field  deployments,  dark  offsets 
were  measured  in  situ  approximately  daily.  Standard  deviations  of  the  change  in  dark  offsets  in 
seawater  within  any  cruise  period  were  significantly  better  than  the  changes  observed  between  air  and 
water.  As  a  percent  relative  to  approximate  count  levels  measured  in  the  lab  for  c=0.1  m"'  AZRD, 
standard  deviations  were  33%  and  35%  worst  case  for  the  Hawaii  cruises.  This  improved  substantially 
for  the  North  Atlantic  cruise  (9%  worst  case)  and  the  Hudson  cruise  (1.7%  worst  case).  The  Hudson 
cruise  was  the  only  period  where  isolated  battery  power  was  used  exclusively,  and  this  presumably 
contributed  to  the  excellent  replicability  in  dark  offsets.  A  conclusion  is  that  working  off  battery 
power  with  a  consistent  sensor  suite  configuration  can  provide  sufficiently  reproducible  dark  offsets. 

The  relationship  between  MASCOT  counts  in  any  channel  and  the  particulate  scattering  at  650  nm, 
bp650,  is  expected  to  follow: 

[COUNTS  -  DO]  =  M  *  bp650  *  exp  (-  L  *  (cp650  *  F  +  W)),  (1) 

where  DO  is  the  dark  offset,  M  is  a  linear  slope  that  is  related  to  relative  gain,  L  is  the  pathlength  (0.2 
m),  cp650  is  particulate  attenuation,  F  is  the  fraction  of  particulate  attenuation  that  is  applicable  (i.e., 
we  do  not  want  to  include  the  scattered  light  that  is  collected  within  the  sample  volume  or  detector 
aperture),  and  W  is  the  pure  water  absorption  (water  scattering  is  assumed  negligible)  at  650  nm  (0.34 
m'^).  Fig,  3  shows  an  example  relationship  between  counts  and  cp650  (assumed  approximately  equal 
to  bp650)  for  AZRD  in  one  of  the  lab  experiments.  As  attenuation  increases,  the  rollover  due  to 
attenuation  along  the  path  is  evident.  Detector  saturation  at  -64000  digital  counts  (14  bit)  is  observed 
in  some  channels. 

In  Eq.  1,  there  are  2  unknowns,  M  and  F.  In  theory,  M  and  F  could  be  solved  for  each  channel  using  a 
least-squares  minimization  of  residuals.  There  are  relatively  few  experimental  points,  though,  so  errors 
can  be  expected  to  be  high.  However,  we  know  that  F  should  be  approximately  constant  for  all  the 
channels  (with  the  exception  that  10,  20,  30,  and  40  degree  channels  have  narrower  detector  apertures). 
Thus,  the  normalized  shape  of  each  one  of  the  channel  regressions  should  be  the  same.  We  therefore 
solved  for  F  by  normalizing  all  of  the  channel  regressions  to  area,  essentially  removing  the  slope  (M) 
or  gain  effect.  Using  a  least-squares  minimization  of  the  model  in  Eq.  1  to  all  the  area-normalized  data 
(included  data  from  every  channel)  resulted  in  a  value  of  0.56  for  E.  Thus,  approximately  44%  of  the 
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attenuation  effect  is  being  removed  by  the  collection  of  secondarily  scattered  light  along  the  path  by 
the  sample  volume  and  detector  aperture. 

If  we  invert  Eq.  1  to  solve  for  bp650,  we  obtain: 

bp650  =  [COUNTS  -  DO]  /  M  *  exp  (  L  *  (cp650  *  F  +  W)),  (2) 

If  the  phase  function  of  a  particle  solution  is  known,  P(0),  we  can  then  obtain  P(0): 

p(0  =  [COUNTS  -  DO]  *  P(0)  /  M  *  exp  (  L  *  (cp650  *  F  +  W)).  (3) 

The  term  [  P(0)  /  M  ]  is  the  scaling  factor  (SF).  For  microspherical  bead  solutions,  values  of  P(0)  can 
be  determined  theoretically  using  Mie  theory  and  values  of  M  can  be  determined  experimentally.  With 
F  in  hand,  we  solved  for  M  individually  for  each  channel  using  a  least-squares  minimization  of  the 
model  in  Eq.  1 .  This  was  carried  out  with  a  series  of  microspherical  bead  solutions  for  which  the 
phase  function  may  be  computed  from  theory. 

The  weighting  functions  of  the  17  MASCOT  detectors  were  computed  analytically  using  code  written 
by  Ron  Zaneveld  and  are  shown  in  Fig,  4.  These  weighting  functions  convolved  with  the  P  calculated 
for  the  beads  are  shown  in  Fig,  5.  Scaling  factors  for  each  channel  were  then  computed  by  dividing 
these  P(0)  values  by  the  experimental  M  values. 

Applying  the  calibration  derived  SF  and  F  parameters,  along  with  measurement-specific  DO’s  and 
cp650,  P(0)  computed  using  Eq.  3  for  the  different  suspensions  from  the  AZRD  experiment  are  shown 
in  Fig,  6.  ECO  VSF  measurements  calibrated  during  the  same  bead  experiment  show  very  good 
agreement.  Applying  the  calibration  derived  SF  and  F  parameters,  along  with  measurement-specific 
DO’s  and  cp650,  calibrated  P(0)  obtained  with  the  MASCOT  in  Hawaii  onboard  the  RW  KW  are 
plotted  in  Fig,  7.  Concurrent  ECO  VSF  measurements  are  overlaid  with  circles.  The  ECO  VSF 
measurements  show  excellent  agreement  with  the  MASCOT  data.  This  is  clear  when  looking  at  a 
vertical  profile  from  Hawaii  (Fig,  8).  Preliminary  analyses  with  Fournier-Forand  phase  function  fits 
have  also  shown  that  the  measured  VSFs  may  be  represented  well  by  a  particle  population  with  a 
Jungian  type  particle  size  distribution  with  low  bulk  refractive  index,  although  some  overestimation  at 
mid-angles  is  evident.  It  appears  that  the  combination  of  two  or  more  such  populations  may  be 
required  to  achieve  a  precise  fit,  not  an  unexpected  finding. 

IMPACT/APPLICATIONS 

Progress  and  results  represent  important  steps  toward  the  development  of  a  multi-angle,  in- water  VSF 
device.  Knowledge  of  the  Inherent  Optical  Properties  including  the  VSF  can  be  used  to  predict  and 
optimize  the  performance  of  a  host  of  Naval  operations  that  rely  on  divers,  cameras,  laser  imaging 
systems,  and  active  and  passive  remote  sensing  systems.  These  include  mine  countermeasures,  harbor 
security  operations,  debris  field  mapping,  anti-submarine  warfare,  and  search  and  salvage  operations. 
These  measurements  may  also  be  used  in  environmental  monitoring  and  research  applications  for 
determining  particle  composition. 


5 


Figure  2.  MASCOT  VSF  device  mounted  with  other  optical  sensors  in  a  custom  cage 
designed  to  sample  the  subsurface  domain. 
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Figure  3.  Digital  response  of  MASCOT  scattering  channels  as  a  function 
of  c(650)  measured  with  an  ac-9  in  numerous  suspensions  ofAZRD. 
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Figure  4.  Analytically  computed  weighting  functions  for 
the  1 7  MASCOT  scattering  measurements. 


Figure  5.  Phase  function,  P(  0),  of 1. 992±0. 025  um  beads.  Circles  are 
values  obtained  from  convolving  the  phase  function  with  the  weighting 
functions  for  the  MASCOT  scattering  measurements. 
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Figure  6.  Fully  calibrated  MASCOT  VSFs  obtained  in  suspensions  of  AZRD.  Concurrently 
collected  ECOVSF  VSF  data  are  plotted  using  (+)  symbols. 


Concurrently  collected  ECOVSF  data  are  plotted  as  circles. 
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Figure  8.  (A)  A  profile  from  Hawaii  of  the  17  P(6)  measured  with  the 
MASCOT;  (B)  selected  P(100),  P(120),  and  P(150)  measured  with  the  MASCOT 
and  overlaid  /3(100),  P(125),  and  P(150)  concurrently  measured  with  an  ECOVSF. 

The  ‘M’  and  ‘E’  labels  in  the  legend  refer  to  MASCOT  and  ECOVSF,  respectively. 
Measurements  at  120  and  150  degrees  are  offset  by  0.0002  and  0.0004  m'^ sr\  respectively, 
for  clarity.  Agreement  between  the  MASCOT  and  ECOVSF  measurements  is  excellent. 
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TRANSITIONS 


We  expect  that  our  efforts  in  developing  an  in-water  VSF  device  and  associated  inversion  techniques 
to  better  understand  particle  dynamics  in  natural  waters  will  lead  to  transition  as  operational  tools  for 
the  fleet  and  the  oceanographic  research  community  in  the  future. 

RELATED  PROJECTS 

This  effort  is  related  to  several  ongoing  efforts  to  develop  optical  sensors  and  associated 
biogeochemical  inversion  techniques.  Related  projects  include; 

•  resolving  the  optics  and  dynamics  of  subsurface  bubble  populations  in  the  Southern  Ocean, 

•  developing  novel  harbor  security  monitoring  capabilities  with  Chuck  Trees  and  Jim  Mueller, 

•  developing  improved  vicarious  calibration  and  validation  methods  for  ocean  color  satellite  remote 
sensing, 

•  investigating  the  sources  of  backscattering  in  natural  waters, 

•  developing  tools  for  ocean  observing  systems,  and 

•  developing  a  surfzone  optical  mine  countermeasure  drifter. 
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